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Abstract 

The semi-relativistic equation is cast into a second-order Schrodinger-like equation 
with the inclusion of relativistic corrections up to order (w/c)^. The resulting equation 
is solved via the Shifted-Z expansion technique, which has been recently developed to 
get eigenvalues of relativistic and non-relativistic wave equations. The Coulomb, Os- 
cillator, and the Coulomb-plus-linear potentials used in qq phenomenology are tested. 
The method gives quite accurate results over a wide range of r and any choice of 
quantum numbers n and I. However, a comparison of the present work with those of 
Lucha et al. and Nickisch et al. will serve as a test of this approach. 
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1 Introduction 



The problem of solving the semi-relativistic equation (SR) is a very old subject, but it 
still an important theme in the existing literature. In recent years much attention has been 
focused on the semi-relativistic equation, because it help in understanding the structure 
of heavy- light mesons, especially the decay constant of the B meson Jb, and the Cabibbo- 
Kobayaski-Maskawa (CKM) matrix elements Vcb and Vub- 

This spatial wave equation (SR), is the more realistic starting point for the descrip- 
tion of bound states consisting of scalar bosons as well as of the spin-averaged spectra of 
bound states consisting of two fermions of masses mi, m2 and total energy M. In the 
past few years, the (SR) equation has been tested for its bound-state energies numerically 
and analytically using different techniques [2] and references therein. All these attempts to 
obtain the energy eigenvalues and corresponding wave functions of bound states from (SR) 
equation rely very heavily on numerical methods. Here, this work formulates the Shifted-/ 
expansion technique (SLET) in which such awkward difficulties may not arise. Very re- 
cently we introduce the Shifted-/ expansion technique (SLET) to solve the Schrodinger 
[3], Dirac [4], and Klein-Gordon equations for some model potentials [5]. SLET simply 
consists of using l/l as an expansion parameter, where I = I — fi. / is the angular mo- 
mentum, and ^ is a suitable shift introduced partly to avoid the trivial case / = 0. In 
addition to that, ^ is chosen so that the next leading term in the energy eigenvalue series 
vanishes, and therefore improving the convergency of the perturbation series. This choice 
is physically motivated by requiring agreement between SLET analytical expressions for 
the eigenvalues and the modified eigenfunctions to reproduce the exact analytical expres- 
sions of the harmonic Oscillator and Coulomb potentials [6,7]. Unlike other perturbation 
methods [8-11], SLET puts no constraints on the coupling constants of the potential or on 
the quantum numbers involved. Therefore, the main motivation of the present work is to 
overcome the shortcomings of previous approaches and to formulate an elegant algebraic 
approach to yield a fairly simple analytic formula which will give a rapidly converging 
leading-orders of the energy values with good accuracy. It is worth mentioning that SLET 
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not only reproduces sufficiently accurate energy values for coulomb plus linear potential 
but it also gives exact analytic results for the ground- state energy coulomb problem. 

In this spirit, this paper is organized as follows. Section 2, recalls the conversion of 
th SR equation into an equivalent Schrodinger- type equation, and there the Shifted-/ 
expansion technique (SLET) for this equation with any spherically symmetric potential is 
introduced. The analytical expressions for SLET's are cast in such a way that allows the 
reader to use them without proceeding into their derivation. Finally, Sec. 3. being devoted 
to test our technique for Coulomb, Oscillator interactions, and for the realistic Coulomb- 
plus-linear potential used in qq phenomenology, and then we conclude and remark therein. 

2 SLET for the Semi-Relativistic Equation with Spherically 
Symmetric Potentials. 

The (SR) equation is the combination of relativistic kinematics with some static 
interaction potential [1], 



where Eni = M — mi — m2, jj, = mim2/{mi + 7712) is the reduced mass and v = 
mf m2/(m^ + 7712) is a useful parameter. The first two terms and (7711+7712) follow from the 
expansion of (p^ +77if )^/^ + (p^ +771! y^"^. In this reduction formalism the spherically sym- 
metric potential V{r) which represents the interaction between the two particles remains 
unspecified. In the physical literature, the second term in Eq.(2) is treated as perturba- 
tion by using trial wave functions [12]. In this work, in order to obtain a Schrodinger- like 
equation, this term is treated using the reduced Schrodinger equation [13]: 



[(p2 + mi fl^ + ip" + ml ) V2 + v{r) - M ] ^(r) = . 



(1) 



An expansion in the powers of {v/c)'^ up to two terms yields 




(2) 



= 4|U^ M - mi - m2 - V{r) 



(3) 
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Thus, the second order Schrodinger— hke equation to order {v/c)'^ becomes 
where r] = I'/fJ-^- Writing the operator in the spherical polar coordinates 



P 



2 ' 



(4) 



(5) 



c?r^ r dr 

and for states of definite orbital angular momentum define the reduced radial wave 
function Rni{r) by ipni{r) = Rne{r)Yim{d , <t>)-, and after some manipulation Eq.(4) can 
be transformed to a Schrodinger- like equation ( in units /i = c = 1) for the radial wave 
function Rni{r) 



1 d2 1(1 + 1) Er,iV{r) 



Rndr) = + Ene)Rne{r), (6) 
2r] 



with 7(r) = V{r) - V'^{r)/2r). 

If the angular momentum / is shifted through the relation I = 1 + Eq.(6) becomes 



^_cP_ ^ [P + 1{2P + 1) + PiP + 1)] ^ ^ ^»|V^(r) 



2nr^ 



F2 



(7) 



where n in this paper is the radial quantum number. 

The systematic procedure of SLET begins with shifting the origin of the coordinate 
through the definition 

2; = /i/2(r-r„)/r„ (8) 

where r„ is currently an arbitrary point to perform Taylor expansions about, with its 
particular value to be determined below. Expansions about this point yield: 
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j=0 



dri 



(9) 
(10) 
(11) 



It should be mentioned here that the re-scaled coordinate, Eq.(8), has no effect on the 
energy eigenvalues, which are coordinate-independent. It just facilitates the calculations 
of both the energy eigenvalues and wave functions. It is also convenient to expand E^t as: 



Enl — ^ Eil \ 



(12) 



2 = 



Substituting Eqs.(9)-(12) into Eq.(7), and expanding around a; = in powers of x and 
one gets; 



2fj, dx^ 2fj, 



I 



f X , Yiro^oX i'{ro)rlx'^ i"{ro)rlx^ 



21 



+ 



+ ■■■) 



vQ 



(l>nl{x) 



(13) 



where the primes of V{ro), and j(ro) denotes derivatives with respect to r^. 



o 

Q 



2ri T] 



EqEi 



E0E2 El 1 
r] 2r]'l 



/ „ EnE-i E-\ Eo . 1 
r] r] F 



(14) 



and Q is a constant that scales the potential, and is set, for any specific choice of T, equal 
to P at the end of the calculations. This is because, the angular momentum barrier term 
goes like P at large T, so should the potential V{r). Equation (13) is exactly of the type 
of Schrodinger-like equation for the one-dimensional anharmonic oscillator and has been 
investigated in detail for spherically symmetric potentials by Imbo et al [6]: 



(15) 
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To estimate P{x) and eg. one has to match the terms of Eqs. (13-14) with the terms in 
Eq.(15) consequently, the final analytic expression for the energy eigenvalues appropriate 
to (SR) particle is 



I 



1 ^ riV{r„)Eo ^ rg7(r,) 



2n r]Q 



Q 



+ 



+- 



H 7^ h Qfi 



2ju 



+ 



rlV{ro)E^ 



+ Q!2 



(16) 



where ai and a2 appearing as a correction to the leading order of the energy expression 
are defined and listed in the appendix of Ref. [4], and the relevant quantities e's and 6' s 
corresponding to (SR) equation are given in the appendix of this paper. A quantitative 
estimate for the energy terms in Eq.(16) can be obtained by comparing the terms of same 
order in / of Eq.(14). The calculations, and the results are very simple. 



V M'^o 

Here, /3 is chosen so that the next contribution to the leading term in the energy eigenvalue 
series to vanish, i.e., Ei = 0, which implies that 

P = -1/2 - fj.{n + l/2)uj, (20) 
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where 

U) = 

and Q satisfies 



a; = - [3 + roF"(ro)/F'(ro) - fxr'oV (rof / (Qr^)] , (21) 



Q = ^[riV'{ro)]\l + 0, (22) 



with 



e = ^1 + [27?/roF'(ro)]^ (23) 



On the other hand, ro is chosen to minimize £^0, such that, 



— ^ = and — ^ > , 24 
dro dr^ 



collecting these terms, and carrying out the mathematics, one can get 

l + 2i + f,{2n + l)ij = r^F'(ro) + MilV^', (25) 



which is an explicit equation in ro- Substituting the expressions of Eq.(17), Eq.(18) and 
Eq.(19) into Eq.(12) immediately gives an expression for the energy eigenvalues, that is 



0^(1) Q;(2) 



(26) 



It is convenient to summarize the above procedure in the following steps: (a) Calculate ro 
from Eq.(25) and substitute it in Eq.(22) to find Q (b) Substitute Q in Eq.(17) to obtain 
Eg. (c) Finally, one can obtain Eg and then calculate Eni from Eq.(26). However, one 
is not always able to calculate ro in terms of the potential coupling constants since the 
analytical expressions become algebraically complicated, although straightforward. There- 
fore, one has to appeal to numerical computations to find rg and hence Eg. 
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3 Applications, results, and discussion. 



In the existing literature, considerable interest has been attached to the analytic solution 
of the semi-relativistic Coulomb ground-state energy problem. Such ground-state energy 
have been derived in Ref.[14,15], and it is given as 

2 m 



E„ 



n = 0, 1,2, 



(27) 



[1 + 12;^] 

On the other hand, an analytical upper bound on the ground-state energy level have been 
obtained by some straightforward variational procedure [16], which reads 



E„ < 2m 



1 



a 



1/2 



n 



0,1,2,... 



(28) 



(2n + 2)2 

In this work, in order to remove this inadequacy, we propose the Shifted/ expansion tech- 
nique, SLET, to solve the (SR) equation, which not only restores the exact analytic expres- 
sion for Coulomb ground-state energy problem, but also predicts quite accurate results for 
the Oscillator, and Coulomb-plus-linear potentials. Therefore, it is worthwhile to illustrate 
that our formulas Eq.(17), (20), (21), and (22) provide a remarkable accurate and simple 
analytic expressions to the S-wave (SR) equation energy eigenvalues of spherically sym- 
metric Coulomb potential V{r) = —a/r. This test is stringent to this approach. Following 
the SLET procedure one obtains uj = 2/m, with m = mi = m2, 

Q = (2n + 2)2/4, 



and 



(2ra + 2f 
2m 



1/2 



(29) 



(30) 



.q;2 (2n + 2)2. 

Substitution of the above expressions into Eq.(17), immediately gives an expression for 
the S-wave eigenvalues 



Eo = 2m 



a 



hence, 



M = 2m 



a 



(2n + 2)2 
1/2 



1/2 



(2n + 2)^ 



n 



2 m. 



0,1,2,. 



(31) 



(32) 
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In Table 1, where the / = eigenvalues obtained from SLET are compared with the 
coulomb analytical result obtained by [14,15], and by some straightforward variational 
procedure [16]. It is found that, the eigenvalues obtained from Eq.(31) are in good agree- 
ment with the results obtained by [14,15], and match very well an upper bound derived 
by [16]. 

In a second test, Eq.(26) is solved as stated in Sec. 2, it is shown that the solutions 
for the Oscillator potential are very well behaved. In Table 2, the convergence of the 
eigenvalues with increasing n, and with the results obtained in Ref.[2] are quit well. 

Finally, the rate of convergence of the eigenvalues obtained from SLET for a realistic 
Coulomb-j^lus-linear potential are shown in Table 3 compared to other methods. 

As a result, it is observed that, the convergence of the eigenvalues are very rapid as 
n, and / are increasing. Therefore, we can say, this technique is easy to implement, the 
results are sufficiently accurate for practical purposes, and the SLET equations are cast in 
away to be used to the unequal mass cases qQ. Moreover, substantial computation time 
reduction has been achieved. 
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Table 1: Energy eigenvalues (in units of GeV) of the {SR) equation, Eq.(31), for the 
Coulomb potential V{r) = —o./r with a = 0.25, and mi = m2 = m = 1.45 GeV. The 
values in round parentheses are those of square root method, whereas, the values in square 
parentheses are those of integral method Ref. [2]. 



n=0 



n=5 



En 
En 
En 

En 



-0.022394 
-0.02274 
1-0.02306) 
[-0.02251] 



-0.005648 
-0.00566 
(-0.00574) 
[-0.00556] 



-0.002514 
-0.002516 
(-0.00254) 
[-0.00244] 



-0.001415 
-0.001415 
^-0.00143) 
[-0.00140] 



-0.000906 
-0.000906 
(-0.000912) 

[-0.00085] 



-0.000629 
-0.000629 
(-0.000635) 
[-0.00065] 
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Table 2: Energy eigenvalues (in units of GeV) of the (SR) equation, Eq.(26), for the 
Oscillator potential V{r) = with mi = m2 = m = 1.310 GeV. The values in round 
parentheses are those of square root method [2], whereas, the values in square parentheses 
are those of Miller's method Ref. [2, 17]. 



(SLET) 





n=0 


n=l 


n=2 


n=3 


n=4 


£=0 


1.6536 


3.5048 


5.1409 


6.6269 


8.0049 




(1.6595) 


(3.5280) 


(5.1654) 


(6.6553) 


(8.0395) 


Miller 


[1.6595] 


[3.5280] 


[5.1654] 


[6.6553] 


[8.0394] 


i=l 


2.6609 


4.3719 


5.9218 


7.3484 


8.6823 




(2.6663) 


(4.3932) 


(5.9441) 


(7.3737) 


(8.7125) 


Miller 


[2.6663] 


[4.3932] 


[5.9441] 


[7.3737] 


[8.7124] 


£=2 


3.6066 


5.2086 


6.6844 


8.0577 


9.3508 




(3.6110) 


(5.2268) 


(6.7039) 


(8.0796) 


(9.3766) 


Miller 


[3.6110] 


[5.2268] 


[6.7039] 


[8.0796] 


[9.3765] 



Table 3: Energy eigenvalues (in units of GeV) of the {SR) equation, Eq.(26), for V{r) = 
-a/r + (3r with a = 0.25, /3 = 0.18 GeV"^, and nii = m2 = m = 1.45 GeV. The values in 
round parentheses are those of square root method Ref. [2]. 





n=0 


n=l 


E^i (SLET) 
n=2 


n=3 


n=4 


£=0 


0.4930 


1.0069 


1.3988 


1.7323 


2.0295 




(0.4924) 


(1.0022) 


(1.3925) 


(1.7252) 


(2.0205) 


£=\ 


0.8342 


1.2484 


1.5971 


1.9053 


2.1855 




(0.8345) 


(1.2481) 


(1.5960) 


(1.9033) 


(2.1793) 


£=2 


1.0958 


1.4600 


1.7796 


2.0685 


2.3345 




(1.0962) 


(1.4601) 


(1.7797) 


(2.0687) 


(2.3346) 
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Appendix 

The definition of e'jS and 6[s for ai and 02 in Ref.[4] are given here which are applicable 
to the {SR.) equation 

.1 = ^^, . = 1,2,3,4. ^. = ^^, . = 1,2,3,4,5,6, (33) 



where 



£i = -(2^ + l)^ , £2 = 3(2^ + l)/2^, (34) 



£3 = -2/m + ^ [7"' (ro) + V" {ro)Ejr,] , (35) 



6Q 



£4 = 5/2^/ + ^ [7"" (rj + F"" {ro)Eo/rj] , (36) 



^i = -^(^ + 1)/m+ ^°^^'^°^^^ (37) 



62 = 3/3(/3 + 1)/2m + ° ^ (38) 



^3 = -2(2^ + 1)/^ , ^4 = 5(2/3 + l)/2^, (39) 



^5 = -3/m + [7 {To) + V {ro)Eo/v] , (40) 
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g 

r ////// y ^ ////// , ^ , 1 y 

^6 = 7/2^ + ^(7 {ro)+V {ro)Eo/v]. (41) 
The terms including Ei have been dropped from the above expressions since Ei=0. 
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